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We present results for the light quark masses for the Wilson quark action obtained with the PCAC relation for 
the one-link extended axial vector current in quenched QCD at f3 — 5.9 — 6.5. This method leads to a remarkable 
improvement of scaling behavior of the light quark masses compared to the conventional method. We obtain 
rni — 3.87(37)MeV for the averaged up and down quark mass and ms ~ 97(9)MeV for the strange quark mass in 
the MS scheme at = 2GeV. 



1. Introduction 

The chiral symmetry breaking term in the Wil- 
son quark action causes large scaling violation ef- 
fects of 0{a) for physical quantities in numeri- 
cal simulations of lattice QCD. A manifestation 
of this effect is that the light quark masses for 
the Wilson action show a strong a dependence at 
finite lattice spacings in contrast to a weak de- 
pendence for the case with the Kogut-Susskind 
(KS) action which retains U(l) chiral symme- 
try [Q. Uncertainties in a long extrapolation to 
the continuum limit needed for the Wilson case 
forms a part of the difficulty to settle the question 
whether the Wilson and KS quark actions give a 
consistent result in the continuum limit This 
leads us to reconsider the definition of quark mass 
for the Wilson action. 

It is well known that the current quark mass 
defined by the PCAC relation for the Wilson 



*presented by Y. Kuramashi 



action Ig,^ differs by 0{a) from the convention- 
ally defined quark mass at finite lattice spacing. 
We examine this point in more detail and find 
that the former definition applied with the one- 
link separated axial vector current that naturally 
arises in chiral Ward identities |^ leads to a sig- 
nificantly improved scaling behavior for the light 
quark masses. 

Our calculations are carried out as a part of our 
simulation for the kaon i?-parameterQ with the 
Wilson quark action in quenched QCD. The pla- 
quette gauge action is employed, and data are col- 
lected at (3 = 5.9 — 6.5. Point source quark prop- 
agators calculated for four values of the hopping 
parameter at each (3 are used for this work. Errors 
are estimated by the single elimination jackknife 
procedure for all measured quantities. 

2. Calculational method 

The conventional(CV) definition of quark mass 
is given by m^^a = {l/K - l/K^^)/2 where 
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Figure 1. rn^ and p as a function of 1/{2K) at 
(3 = 6.3. 

the critical hopping parameter K^^ is deter- 
mined by extrapolating the pseudo scalar me- 
son mass squared mjjg linearly to mj^g = 
in the inverse hopping parameter 1/{2K). To 
extract the averaged up and down quark mass 
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)/2 and the strange quark 
mass m^^ we assume the following quark mass 
dependence for mps and the vector meson mass 
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The lattice spacing a is fixed with using 
nrip = 770MeV. We employ the hadron mass ra- 
tio mT^jmp = 0.1783 to determine the averaged 
up and down quark mass mp^. The strange 
quark mass is estimated in two ways using 
TOif/wp = 0.644 and m^/mp = 1.323. We con- 
vert rn^ into m^^ defined in the MS scheme at 
the scale ^ = 1/a using the tadpole- improved per- 
turbative mass renormalization factor ZmifJ-a = 
1) evaluated with aj^(l/a). 

The Ward identity (WI) method requires a cal- 
culation of the p parameter defined by the PCAC 
relation 01; 



2p{K) = 



|V^A°"*-°|7r'^(p = 0)) 
(0|P'^|7r«(p = 0)) '■ 



(3) 



We take the one-link extended axial vec- 
tor current A™*^" since it is this current 




Figure 2. Averaged up and down quark mass in 
the MS scheme at /i = 2GeV as a function of a. 

which naturally arises in the Ward identities. 
We extract the p parameter from the ratio 

(V4Af*(t)p(i))/(P(t)P(i))/4+(V4Af*(r-t+ 

2)P(l))/(P(r - t + 2)P(l))/4, each two-point 
function projected to the zero spatial momentum, 
by fitting a plateau as a function of t. 

In Fig. a representative result for the p pa- 
rameter is plotted as a function of 1/{2K) to- 
gether with m^. We observe a clear linear be- 
havior both for p and mj. The critical hopping 
parameter extracted from a linear extrapo- 
lation for p is slightly different from K^^ for 
(see Table |l] for numerical details). We ascribe the 
discrepancy to uncertainties in the extrapolations 
of p and , because the critical hopping parame- 
ter obtained with the two definitions should agree 
at each /3. To avoid this problem, we define the 
bare quark mass for the WI method m^^{l/a) by 
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(1/a) 



(4) 



We convert quark masses defined on the lattice 
into those in the MS scheme at the scale p — 1/a 
using the tadpole-improved perturbative renor- 
malization factor Zji<sxt[pa — 1)/Zp(pa — 1) 
evaluated with a-^{l/a). 

We quote final results for m^^'^^ at the scale 
p — 2GeV which are obtained by a two-loop 
renormalization group running from p = 1/aQ. 
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Table 1 



Critical hopping parameters and the light quark masses in the MS scheme at ^ — 2GeV. 







input 


/3=5.9 


6.1 


6.3 


6.5 


a = 








0.15986(3) 


0.15502(2) 


0.15182(2) 


0.14946(3) 










0.15976(1) 


0.15497(1) 


0.15179(1) 


0.14948(1) 






(MeV) 


m^/mp 


5.54(14) 


5.14(22) 


4.89(30) 


4.62(41) 


3.96(49) 




(MeV) 


m-„/mp 


3.97(10) 


3.94(17) 


3.97(23) 


3.86(32) 


3.87(37) 




(MeV) 


■niK/mp 


139(3) 


129(5) 


123(8) 


116(10) 


100(12) 




(MeV) 


rriK/mp 


100(3) 


99(4) 


100(6) 


97(8) 


97(9) 




(MeV) 


m4,/mp 


184(13) 


149(17) 


151(24) 


120(23) 




— WI 

mr 


(MeV) 


m^/mp 


132(9) 


115(13) 


122(20) 


99(19) 





3. Results for light quark masses 

In Table [| we summarize the results for the av- 
eraged up and down quark mass m;. The lattice 
spacing dependence of m/ is shown in Fig. ^. The 
results with the WI method show a remarkably 
flat behavior. This allows a reliable extrapola- 
tion to the continuum limit by a linear function 
in a, with which we find mf^ = 3.87(37)MeV at 
a = 0. In contrast a large scaling violation effect 
is seen for the results with the CV method. While 
a linear extrapolation yields a consistent result 
= 3.96(49)MeV in the continuum limit, sys- 
tematic uncertainties due to a long extrapolation 
are quite large. Finding a theoretical explana- 
tion as to why finite a corrections are so small for 
quark masses with the WI method would be an 
interesting problem Q. 

In Fig. 1^ we also plot the KS results |^ for com- 
parison. They show small scaling violations and 
are systematically smaller than the results for the 
WI method. However, we cannot conclude at this 
stage whether the Wilson results are inconsistent 
with the KS results. The one-loop correction in 
the mass renormalization factor for the KS quark 
action is 50% — 100% depending on the lattice 
spacing for the results in Fig. ||, which suggests 
that higher order corrections might be large. 

Our results for the strange quark mass TOs de- 
termined from rriK and are listed in Table 0. 
We find that results with the WI method are quite 
flat also in this case. 

Estimates of rris using rriK are expected to 
satisfy Wg ~ 25m; for each (3 if we assume 
the functional form (|l|), since then m\/m'^ = 



{nil + Wis) / {2fni) ~ 13. Making linear extrap- 
olations in a for the results of the WI and the 
CV method we find a mutually consistent result 
in the continuum limit: = 97(9)MeV and 

mjv = 100(12)MeV. 

For the alternative determination of Wis with 
the use of m^, the discrepancy between the WI 
and the CV results reduces toward the contin- 
uum limit. However, large errors originating from 
those of the vector meson mass hinder us from re- 
liably extrapolating the results to the continuum 
limit. 
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